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ABSTRACT 



Context. Chromospheric activity has been calibrated and widely used as an age indicator. However, it has been suggested that the 
viability of this age indicator is, in the best case, limited to stars younger than about 1.5 Gyr. 
Aims. I aim to define the age range for which chromospheric activity is a robust astrophysical clock. 

Methods. I collected literature measurements of the S-index in field stars, which is a measure of the strength of the H and K lines of 
the Ca II and a proxy for chromospheric activity, and exploited the homogeneous database of temperature and age determinations for 
field stars provided by the Geneva-Copenhagen survey of the solar neighbourhood. 

Results. Field data, inclusive data previously used to calibrate chromospheric ages, confirm the result found using open cluster data, 
i.e. there is no decay of chromospheric activity after about 2 Gyr. 

Conclusions. The only existing indication supporting the viability of chromospheric ages older than 2 Gyr is the similarity of chro- 
mospheric activity levels in the components of 35 dwarf binaries. However, even in the most optimistic scenario, uncertainty in age 
determination for field stars and lack of sufficient data in open clusters make any attempt of calibrating an age activity relationship for 
old stars premature. The hypothesis that chromospheric activity follows the Skumanich law, i.e. that it is proportional to t"''^, should 
be relaxed. * 
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1. Introduction 

Chromospheric activity (CA) is the amount of energy emitted in 
the chromosphere, i.e. the layer of a late-type star above the tem- 
perature gradient inversion. This energy is radiated especially 
through emissions of strong resonance lines, whose flux is there- 
fore a proxy for CA. The S index is a measure of the strength of 
the chromospheric emission core of the H and K lines of the 
Ca II, F^^ is its conversion into flux, and R^^^ is the flux nor- 
malized to the bolometric emission of the star. These are the 
most commonly used proxies for CA, and they have been mon- 
itored for decades in nearby stars in the Mount Wilson (see e.g. 
iBaliunas et al.L [19951) and, more recent l y, in s everal other obser- 
vational campaigns (see e.g. lHall et alllBoTh . 

The CA declines during the stellar life time, as does the rota- 
tion, which powers it. Therefore CA is in principle a good astro- 
physical clock and, as such, has been calibrated and widely used. 
However, CA, on top of its long-term decay, undergoes short- 
term variations. The Sun, for instance, has an activity cycle of 
1 1 years with secular variations that caused the Maunder min- 
imum. This sets limits on the precision of chromospheric ages 
based on one-epoch CA measurements. 

Observations in open clusters suggested that CA in stars 
in open clusters at about 1.5 Gyr, such as IC465 1 and 
NGC 3680, already dropped to sol ar levels (iPace & Pasquinil 
l2004tlLv ra & Porto de Mell^ l2005l) . This result was confirmed 
by Zhao et alj (i201 ll) who studied wide binaries. iPace et aP 
(2009) claimed that even in the age range between 0.7 and 1.2 
Gyr any CA evolution must be less important than short-term 
variations, and that a sudden drop of CA occurred on a very 



short time-scale : 0.2 Gyr. These con clusions relied on the age de- 
terminations of ISalaris et al.l (1200 4'). In fact, according to these, 
the ages of the active clusters, which all spanned roughly the 
same interval of CA levels, ranged from 0.7 Gyr (Hyades and 
Praesepe) to 1.2 Gyr (NGC5822), while the ages of NGC3680 
and IC4651, which are as inactive as the Sun, were estimated to 
be 1.4 and 1.7 Gyr, respectively. A more recent set of homoge- 
neous age determinations for these clusters makes NGC5822 0.9 
Gyr (Can-aro et al., 201 ]}, IC4651 1.5 Gyr, and NGC3680 1.75 
Gyr old ( Anthony-Twaro g et al.L 20091) . Therefore, according to 
the CA measurements of iPace et al.l and these new age determi- 
nations, the transition from chromospherically active to inactive 
would occur in 0.6 Gyr, and, most important, the interval of time 
during which CA remains virtually constant is reduced to 0.2 
Gyr. In other words, it is still waiTanted to investigate the viabil- 
ity of CA as an age indicator until 1.5 Gyr 

In contrast to the result of Pace & Pasquinil (l2004l) . 
Lvra & Porto de Mellol (l2005h . and IZhaoetal.l (l201lh . several 
calibrations of CA as an age indicator have been published (see 
e.g. Soderblom et al., 1991; Donahu^ 119931: iLachaume et al.L 
[1999: Mamajek & HiUenbrand. 2008) based on field stars, bi- 
nary systems, or open-cluster members, and these calibrations 
are claimed to be valid at all ages. 

The problem that all of these studies face is the paucity of 
good data for open clusters older than about 1 Gyr, M67 being 
the only such cluster with a rich sample of high-quality CA mea- 
surements made for stars spanning a wi de range of temperatures 
(Giampapa et al., 2006). The sample of iPace & Pasquinil (|2004|) 
and Pace et al.i (i2009i) include, in addition to M67 members, only 



1 



G. Pace: Chromospheric activity as age indicator 



seven stars older than 1 Gyr, two dwarfs in NGC3680 and five in 
IC4651, and they are all hotter than the Sun. 

iMamaiek & Hillenbrand (2008) reported only two more 
open clusters older than 1 Gyr with CA measurements: 
NGC752 (15 stars obser ved) and NGC188 (3 star s observed) 
origin ally observed by iBarrv et aP (Il987b and iBarrv et all 
(1984). Ho weve r, thes e dat a were co nsidered suspicious 
(Soderblom etal], |1991|) and iMamaiek & Hille nbrand (2008J 
did not use them either in their Table 5 or in their Fig. 4, where 
other cluster data are reported. 

Clearly, more open-cluster data are needed. However, since 
the publication of the most recent CA calibration I have knowl- 
edge of, more S-index measurements were published, and new, 
precise, homogeneous determinations of stellar ages and tem- 
peratures for the solar neighbourhood have been made available. 
This allows us to reexamine CA evolution through the study of 
field stars. 



2. The data 

2.1. Field data 

All the available data in the literature reporting S-index 
measurements for field dwarf stars in the same scale as Mount- 
Wilson data were collec t ed. The data include the fo llowing 
sources: lArriagadal (1201 ll) . Jenkins et al.l (201 1 ). Baliunas et al.| 



(Il995h. knutson et alJ (l2oToF 



lGravetal.1 (12 003'). 'S chroder et - _ 
(I2007h. iLoc kwood et all (l2007h 



(2010), 



Grav et al.l 



White et al.l 



Lopez-Santiag o et al. 
al.l (I2009'). Cincunegui et al. 
p iDu ncan et al. '(1991), 
iTinnev et air(l2002lK IHalletal.1 (I2007h 



iHenrv et al.l 



and IWright et al.l 



( 2004 ). [Strassmeier et al. (2000) measured absolute H and K 
emission-line fluxes of more than 1000 late-type stars, but they 
did not report S-index measurements. I transformed their flux 
measurements into line strength (HK,,,.^,^) by reverting their 
equations 2) and 3), and then to an equivalent of the S-index in 
the Mount Wilson scale, by calibrating a relationship between 
a^sirass and the S-index for the 329 stars in their sample in 
common with at least one of the other works mentioned above. 
The spread around a linear relationship between HKsrrass and 
S-index is roughly 10 %, co mpatible with CA cycles ari d errors. 
I also included data from iBuccino & MauasI ( l2008l) . which 
are measurements of the strength of the Mg II h and k lines, 
transformed into the S index through a calibration made with a 
set of 1 17 nearly simultaneous observations of Mg II and Ca II 
fluxes of 21 main-sequence stars. Data from Lockwood et al. 
and Baliunas et al. are taken from the Mount Wilson survey and 
cover the entire activity cycle for each target. I adopted them 
whenever possible. For the other stars, when more than one 
measurement was available in the literature, the adopted S-index 
was taken as the average between the highest and the lowest 
individual measurement recorded, which best represents the 
time-averaged CA level, because typically stars are not observed 
during an entire CA cycle. Unfortunately, some authors did 
not repo rt single measure ments but only indicated the average 
S-index. [Hall et al] (l2007 h reported S-index measurements only 
for part of their total sample, for the rest of the stars only R^^^ 
values are given, and the photometric information they used 
to transform S-index into R^^ is not indicated either This 
problem aff'ects 27 stars of their sample that could have been 
used otherwise. Only one of these, HD701 10, was not measured 
in any of the other aforementioned works. 

Although S-index measurements come from a variety of dif- 
ferent sources, they were all calibrated to match the same scale. 



which is the Mount Wilson scale. The most significant source of 
error is the intrinsic variability of CA. The conclusions presented 
here are independent of the choice of the sources, i.e. they hold 
if we discard one or more of the studies from the compilation, 
provided that sufficient data are analysed. 

To guarantee homogeneity of temperature and age evalua- 
tions, this information was taken from the Geneva-Copenhagen 
survey of the Solar neighbourhood (GCS): a set of determi- 
nations of metallicity, rotation, age, kinematics, and Galactic 
orbits for a complete, magnitude-limited, and kinematically 
un biased sample of 16682 nearby F and G dwarf stars 
(No rdstrom et al.L 12004). Two revised versions of the GCS were 
publi shed by the same grou p based on improved uvby calibra- 
tion dHolmberg et al.L l2007h and o n revision of the Hipparcos 
parallaxes (*Holmberg e t al. . l2009 l). I used the last revision by 
Casagrande et al. (201 1), which adopts the infrared flux method 
(IRFM) to improve the accuracy of temperature estimations. For 
this method, with the most up-to-date calibration, the uncertainty 
in the determinations of the Tg/f is of the order of 30 K, plus 
a systematic error, i.e. uncertainty in the zero point, of 15 K 
(Casagrande et al., 2006, 2010). Accurate information on bina- 
rity is also given in this last version of the GCS, and I used it to 
select single stars. 

By cross-correlating the collection of S-index measurements 
with the GCS data set and discarding multiple systems, I finally 
selected 1744 single dwarves. The data collected about these 
stars are available at the CDS in electronic form. From these, 
the 494 stars with age estimations precise enough for the scope 
of the present work were additionally selected icr„gg <2 Gyr or 
^ <0.3). 

age ' 

S-indices were transformed into R^^ following the proce- 
dure of INovcs et al. (1984), which involves B-V colours. The 
adopted values for B-V wher e obtained from the tem peratures 
by reverting the calibration of lCasagrande et aD(l2010l) . 



2.2. Open cluster data 

I used here lo g R'fff evaluations and B-V co lours from the 
compilation of Mamaiek & Hillenbrand! (12008 ^ for the Ursa 
Major moving group (Collinder 285), Hyades, and M67, while 
I ignored data on younger clusters of the same source, since 
they are not relevant for the scope of t he present work. I 
added data of NGC75 2 and NGC188 from iBarrv et"al1 ([T987b . 
and from iPace et al.l (12009) the data of Praesepe, IC4756, 
NGC5822, IC4651, NGC3680, as well as one Hyades mem- 
ber and four M67 members not i ncluded in the sample of 
IMamaiek & Hillenbrand! iPace et al] used spectroscopic temper- 
atures for their stars, and I adopted them here. The other temper- 
atures were obtained through the calibration of Casa grande et al.l 
(.2010) from B-V dereddened colours given in the respective 

source. 

The S-index measurements of iBarrv et all were made on 
spectra at a lower resolution (~2A FWHM) than the classi- 
cal Mount- Wilson data, but they were converted to match the 
Mount- Wi lson system. However, when I transformed them into 
R^^ using lNoves et al.l (11984!) . I obtaine d results com pletely in- 
consistent with the values published in iBarrv et all with typi- 
cal differences of 0.1 in logarithmic scale, and even larger when 
taking into account the correction suggested bv lSoderblom et ali 
(1991). I believe that Barry et al.'s S-index measurements are 
trustworthy, since their conversion into R^^ gives reasonable re- 
sults, while both Ba rry et al.'s value s of R ^^ x 10^ and their 
correction following ISoderblom et al.l (1 199 it are suspicious. As 
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Fig. 1. log RC,^ versus age diagram. 



for the former estimations, in addition to being inconsistent with 
the S-indices from which they are supposed to be derived, the 
highest log R^^^ value for NGC752 exceeds that for the most ac- 
tive star among Hyades and CoUinder 285 by 0.38 dex. This is 
highly implausible since Hyades and Collinder 285 are younger 
than NGC752 and their sample considered here outnumbers by 
more tha n 3 to 1 that of NGC752. By introducing the cor- 
rection of ISoderblom et all this probably spurious value is de- 
creased by only 0.01 dex in the logarithmic scale, and one ob- 
tains log R^^=-5.41 dex for the least active of the three stars 
of NGC188. This is not quite as suspicious as the case for the 
most active star in NGC752, but still calls for some mistrust. 
In fact, it is 0.28 dex lower than the log R^^ value of the least 
active of the 26 stars in M67. In addition, among the 1744 sin- 
gle stars with CA measurements collected here, only 12 stars 
have log R'f^f^ <-5.41 dex, while the most active NGC188 mem- 
ber, with lSoderblom et aLf s correction, would have a CA above 
the average of field stars. I therefore used R^^ values derived 
from published S-indices for NGC752 and NGCI8 8. Because I 
failed to identify the star numbering system used bv lBarrv et al.L 
and they did not report the stellar coordinates, I could only use 
the B-V colours given in the paper without updating the pho- 
tometry and membership information. The CA measurements 
in NGC752 and NGC188 should be taken w i th a gr ain of salt. 
The S-index measurements from Pace et al.l (l2009l) are based 
on high-signal-to-noise and high-resolution spectra taken with 
UVES at the VLT, but they are performed in a diff'erent way, us- 
ing only the K-line. I used the 17 stars, either M67 or Hyades 
members and the Sun, in common between Pace et al. ( 2009) 
and th e compilation of open clusters by Mamaiek & Hillenbrand 
(l2008l) t o find a lin ear relationship between the S-index mea- 
sured bv lPace et al .' and the Mount-Wilson one, and to transform 
the former into the Mount- Wilson scale. The 17 stars in common 
span large part of the entire CA range of the open cluster sam- 
ple, and their deviation from the the linear regression is of only 
~ 5%, therefore the transformation is dependable. 

3. Results and discussion 

3.1. Field stars 

In Fig. [U the sample of 494 stars is plotted in a log R^^ versus 
age diagram. Typical errors in CA reported in the figure are esti- 
mated by looking at the full CA variations for the stars measured 
twice or more. For about 95% of those stars, full variations are 
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Fig. 2. log R^^ versus dereddened B-V colour for open clusters 
and the Ursa Major moving group (Collinder 285, UMa in the 
legend). 



within 0.4 dex, and within 0.2 dex for stars older than 4 Gyr, 
therefore I adopted 0.2 dex as typical error for stars younger 
than 4 Gyr and 0.1 dex for the others. 1 cr errors in age were 
imposed, as a selection criterion, to be either less than 2 Gyr or 
less than 30% of the age. As a result, the typical horizontal error 
bar is of about 1.5 Gyr The L-shaped distribution of the bulk 
of the data points in Fig. [1] strongly suggests that the decay of 
CA stops completely at a relative early age. This age may ap- 
pear to be about 3 Gyr in Fig. [1] as a result of the errors that 
scatter the data points horizontally, but could be much younger 
in reality. It is not possible to understand how CA evolves in this 
first phase, because the uncertainty in the age is of the same or- 
der of magnitude as the time span over which the decay occurs. 
No matter how data are further selected on the basis of the num- 
ber of CA measurements available and on the temperature, there 
is no indication of a CA-age anticorrelation for older stars. The 
vertical spread of the points relative to old stars is compatible 
with being caused by cyclic variations, which proves that CA 
measurements are reliable also for inactive stars. The presence 
of very inactive stars only at young ages is connected with the 
temperature. Indeed, among the 1744 single stars with CA mea- 
surements collected here, the 1 1 least active are all hotter than 

6300 K, and are therefore young. 

Following Rocha-Pinto & M aciell lll998h . I checked the cor- 
relation between metallicity and log R^;^ : the Pearson correla- 
tion coefficient between them is roughly -0.2. 



3.2. Open clusters 



Fig. |2| shows a log R^^ versus temperature diagram. The 
ages reported in the legend are taken from the fo l lowin g 
sources (in order of adopted priority): ICairaro et al.l (20TT]): 
lAnthonv-Twarog et al.1 (l2009h : ISalaris et al.1 (|2004 : iKing et al] 
(I2003h . 

Data points in Fig.|2|can be broadly divided into two groups, 
one with the bulk of the young-cluster members (age < 1.45 
Gyr, filled symbols), with a slightly increasing trend of log R^^ 
with decreasing temperature, and one with the bulk of old clus- 
ter members (age > 1 .5 Gyr, empty symbols), characterized by 
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a more complicated pattern of log R^^ as a function of tem- 
perature. Only three of the 78 members of the former group are 
less active than log R^^=-4.63 dex, and only three out of the 89 
members of the latter group, all M67 members, are more active 
than log R^^=-4.60 dex. The reason why three M67 stars are at 
a Hyades activity level is probably binarity, as confirmed by the 
fact that one of these stars is Sanders 747, whose spectrum taken 
with UVES at the VLT (observing run 66.D-0457) clearly shows 
that it is a double-lined binary. It follows that no more than ~ 5% 
of the single stars should be expected to have an activity atypical 
for their age. Although the CA levels spanned by the two groups 
almost overlap, the average difference between them is as large 
as Alog R^^ ~ 0.4 dex, i.e. twice the typical error for young 
stars and four times as much that for old stars. 

According to the most recent age determinations, the only 
cluster of the active group significantly older than the Hyades 
studied so far is NGC752, which is on average slightly less active 
than Hyades and Collinder 285, even though all NGC752 mem- 
bers, with the exception of one outlier, have CA levels within 
the range encompassed by the bulk of Hyades stars. The drop of 
the average value of log R^^ between the age of NGC752 and 
that of the youngest inactive cluster is ~ 0.3 dex or more. The 
young est inactive cluster, according to lAnthonv-Twarog et alj 
d200^ . is IC4651, only 50 Myr older than N GC752. Theref ore, 
even after the revision of open-cluster ages bv'Carra ro et al.l and 
[Anthonv-Twarog et al., we still have an indication of a very rapid 
CA drop. However, the dependability of the data on NGC752 is 
questionable (see Sect. 12.21 ). 

The dependability of the NGC188 data is also questionable, 
and we have only three stars in a narrow temperature range. It 
is interesting, however, that it is the oldest cluster, by far more 
active than IC465 1 , NGC3680, and M67, thus definitely pointing 
to the absence of CA decay after 1 .5 Gyr 



3.3. Accuracy of the measurements. 

Data shown in Fig. [1] and Fig. |2] are both characterized by the 
lack of old and active stars, as expected. The presence in Fig. 
[1] of young and inactive stars can be explained by the errors in 
age. In Fig. [1] the age limit after which no more active stars 
are present is 3 Gyr, while in Fig. |2]this age limit is 1.5 Gyr 
A gain, this is compa t ible w ith the 1 cr errors in the age quoted 
by ICasagrande et al.l (1201 ll) . which, for the upper-left bunch of 
data points of the diagram in Fig. [T] peak around 0.2 Gyr but 
extend all over to the imposed upper limit of 2 Gyr There is 
a significant diff'erence in the CA limit that sets the distinction 
between active and inactive stars: for open clusters this limit is 
at log R^^=-4.6, while looking at the field data, the same limit 
should be the solar CA level. This last circumstance seems to be 
incompatible with the Sun being in the middle of the CA dis- 
tribution of M67 members. This raises the question whether the 
OC data are on a scale dif ferent than t he fiel d data, but this would 
be surprising, since both lPace et al.l (12009) and lGiampapa et all 
t2006,) found that most of M67 members have C A levels that fall 
in the range spanned by the Sun, they both used high-resolution 
and good signal-to-noise spectra of the Ca II H (only Giampapa 
et al.) and K lines, and they have temperature estimations that 
perfectly match for three solar-type stars, but differ systemati- 
cally by about 100 K for the other (hotter) stars. 

4. Summary and conclusions 

I collected a very large sample of published S-index measure- 
ments that were calibrated to match the Mount- Wilson scale. For 



field stars, I used information on temperatures and ages of the 
very exhaustive and homogeneous determinations of the GCS, 
adopting the last revision by Casagrande et al. (201 1). 

Not only open-cluster data, but also field-st a r dat a 
strongly favour the conclus ion of Pace & Pasauinil (|2004|) . 
iLvra & Porto de Mellol (l2005h . and IZhaoet al.l (1201 lb : there is 
no evolution of chromospheric activity after ~ 1.5 Gyr 

Mamaiek& Hillenbrand (2008) collected data on 35 dwarf 
binaries. The CA levels of their components are strongly corre- 
lated, even if we consider only the systems more inactive than 
log RJy^=-4.60. This is the only piece of evidence supporting a 
correlation between CA and age holding for old stars. The search 
for an alternative explanation for this fact, invoking processes 
typical of binary systems, is definitely warranted, but the possi- 
bility that the CA is a tight, strictly monotonic function of time, 
even though extremely unlikely in the light of the analysis pre- 
sented here, cannot be completely ruled out. According to the 
present analysis, this fu nction, if it exists, i s not like oc r^^^, as 
tentatively suggested b v ISkumanichI (1 1 9721) on the basis of very 
scanty data and widely used ever since its pubblication, and our 
present knowledge of it is too poor for it to be an astrophysical 
clock of any use. 

At present, the only prediction on age we can infer from 
CA is whether a star is younger or older than the open clus- 
ter NGC752, which is estimated to be 1.45 Gyr old, with the 
caveat that there exist young and very inactive stars hotter 
than 6300 K. The situation is likely to improve rapidly due to 
the accumulation of asteroseismic data. Indeed, asteroseismol- 
ogy ha s proved to allow age determination with a good pre- 
cision ('Christense n-Dalsgaard et al.L 1201 Ot iMathur et all 120121; 
Metcalfe et al., 201^^ 
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